Atmospheric 
Chemistry 
and Physics 


Atmos. Chem. Phys., 12, 1-12, 2012 
www. atmos-chem-phys .net/ 12/1/2012/ 
doi:10.5194/acp-12-l-2012 
© Author(s) 2012. CC Attribution 3.0 License. 




Toward a combined SAGE II-HALOE aerosol climatology: 
an evaluation of HALOE version 19 stratospheric aerosol 
extinction coefficient observations 

L. W. Thomason 

NASA Langley Research Center, Mail Stop 475, Hampton, VA 23681, USA 

Correspondence to: L. W. Thomason (l.w.thomason@nasa.gov) 

Received: 15 March 2012 - Published in Atmos. Chem. Phys. Discuss.: 5 June 2012 
Revised: 13 August 2012 - Accepted: 16 August 2012 - Published: 


Abstract. Herein, the Halogen Occultation Experiment 
(HALOE) aerosol extinction coefficient data is evaluated in 
the low aerosol loading period after 1996 as the first neces- 
sary step in a process that will eventually allow the produc- 
tion of a combined HALOE/SAGE II (Stratospheric Aerosol 
and Gas Experiment) aerosol climatology of derived aerosol 
products including surface area density. Based on these anal- 
yses, it is demonstrated that HALOE’s 3.46 pm is of good 
quality above 19 km and suitable for scientific applications 
above that altitude. However, it is increasingly suspect at 
lower altitudes and should not be used below 1 7 km under 
any circumstances after 1996. The 3.40 pm is biased by about 
10% throughout the lower stratosphere due to the failure to 
clear NO2 but otherwise appears to be a high quality product 
down to 15 km. The 2.45 and 5.26 pm aerosol extinction co- 
efficient measurements are clearly biased and should not be 
used for scientific applications after the most intense parts 
of the Pinatubo period. Many of the issues in the aerosol 
data appear to be related to either the failure to clear some 
interfering gas species or doing so poorly. For instance, it 
is clear that the 3.40 pm aerosol extinction coefficient mea- 
surements can be improved through the inclusion of an NO2 
correction and could, in fact, end up as the highest quality 
overall HALOE aerosol extinction coefficient measurement. 
It also appears that the 2.45 and 5.26 pm channels may be 
improved by updating the Upper Atmosphere Pilot Database 
which is used as a resource for the removal of gas species 
otherwise not available from direct HALOE measurements. 
Finally, a simple model to demonstrate the promise of mixed 
visible/infrared aerosol extinction coefficient ensembles for 
the retrieval of bulk aerosol properties demonstrates that a 


combined HALOE/SAGE II aerosol climatology is feasible 
and may represent a substantial improvement over indepen- 
dently derived data sets. 


1 Introduction 

In the stratosphere, where aerosol composition is predom- 
inately mixtures of H2SO4 and H2O, aerosol extinction 
at Halogen Occultation Experiment (HALOE) wavelengths 
is dominated by absorption. Variation with wavelength is 
mostly driven by changes in the imaginary index of refrac- 
tion that are themselves modulated by temperature and rel- 
ative humidity. The extinction measurements generally ex- 
hibit a second order dependence on aerosol size. As such, in- 
frared aerosol extinction coefficient measurements are nearly 
linearly related to total aerosol volume but provide lim- 
ited information on aerosol size distribution. This stands 
in contrast to similar measurements in the visible and near 
infrared where absorption by sulfate aerosol is effectively 
zero and extinction is dominated by the positive wing of 
the aerosol size distribution with limited dependence on the 
smallest aerosol present. This can be seen in Fig. 1 which 
shows Mie extinction kernels scaled to per unit volume of 
aerosol and for measurement wavelength for 75 % H2SO4- 
25% H2O aerosol at 220 K at HALOE and some SAGE 
II/III (the Stratospheric Aerosol and Gas Experiment) mea- 
surement wavelengths (a unitless quantity). A substantial 
history of visible/near infrared extinction coefficient mea- 
surements has been produced by the SAGE-series of instru- 
ments (listed in Table 1). The information content of infrared 
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Fig. 1. This figure shows the Mie aerosol extinction coefficient ker- 
nels (per unit volume of aerosol and wavelength-weighted, a unit- 
less quantity) weighted by wavelength for a 75 % H 2 SO 4/25 % H 2 O 
aerosol at stratospheric temperatures for SAGE II/III and HALOE 
measurement wavelengths. Dotted lines are kernels for SAGE II/III 
wavelengths and solid lines are for HALOE wavelengths. 


Table 1. The Stratospheric Aerosol and Gas Experiment (SAGE) 
series of instruments including the Stratospheric Aerosol Measure- 
ment (SAM II). 


Instrument 

Dates of 
Operation 

Aerosol Extinction Coefficient 
Measurement Wavelengths 

SAM II 

1978-1993 

lOOOnm 

SAGE 

1979-1981 

450, 1 000 nm 

SAGE II 

1984-2005 

386, 452, 525, 1020 nm 

SAGE III 

2002-2006 

386, 448, 521,602, 676, 755, 
868 , 1019, 1545 nm 


measurements (from HALOE) and visible/near infrared mea- 
surements (from SAGE II/III) is sufficiently different that 
the combination of these measurements may substantially re- 
duce the uncertainty of indirectly derived aerosol products 
such as aerosol surface area density (SAD) (e.g., Thomason 
et ah, 2008). 

HALOE is a gas-filter correlation radiometer that uses so- 
lar occultation to measure vertical profiles of a number of 
important trace gas species and multiwavelength aerosol ex- 
tinction coefficient from the upper troposphere as high as the 
thermosphere. It was deployed as a part of the Upper At- 
mosphere Research Satellite from the space shuttle Discov- 
ery on 15 September 1991 and operated through November 
2005 when the mission was terminated. HALOE measure- 
ment species include O3, H2O, CH4, NO, NO2, HC1, HF 
and temperature (via CO2 absorption). Data from this mis- 


Table 2. HALOE aerosol measurement locations with primary gas 
species and other absorbers, a indicates that the species is cleared 
using HALOE derived products (clearing is inferred by the inclu- 
sion of their uncertainties in the aerosol uncertainty as noted in 
Hervig et al., 1995). The symbol b indicates that this species is re- 
moved using climatological values (M. Hervig, personal communi- 
cation, 2012). The symbol c denotes that, while NO 2 is measured by 
HALOE, it is not removed from this aerosol product (E. Remsberg, 
personal communication, 2012 ). 


Wavelength 

(nm) 

Primary 
Gas Species 

Other Absorbers 

2.45 

HF a 

H 2 O a , CH 4 a , C0 2 b , O 3 , N 2 0 

3.40 

HCl a 

CH 4 a , 0 3 a , H 2 O a , N0 2 c 

3.46 

CH 4 a 

N0 2 a , H 2 0\ HCl a 

5.26 

NO a 

N 2 0, C0 2 b , H 2 O a , 0 3 a 


sion were instrumental in the discovery of the water vapor 
tape recorder (Mote et al., 1996) and trends in HC1 (An- 
derson et al., 2000; Jones et al., 2011). In addition to the 
gas species, aerosol extinction coefficient profiles for the up- 
per troposphere through the stratosphere are reported at four 
wavelengths (2.45, 3.40, 3.46, and 5.26 pm). Broadly, these 
data show the immediate aftermath of the June 1991 Pinatubo 
eruption and the long recovery of stratospheric aerosol lev- 
els throughout the 1990s. This is followed by a relatively low 
aerosol loading period between 2000 and the end of the mis- 
sion where only a few minor volcanic events disturb an oth- 
erwise quiescent stratosphere. 

HALOE aerosol extinction coefficient profiles are derived 
as a residual from the gas species retrievals at each wave- 
length using a methodology that is described in Hervig et 
al. (1993) with some additional detail reported in Hervig 
et al. (1996). The aerosol extinction coefficient data is cor- 
rected for the gas species measured at those wavelengths as 
denoted in Table 2. In addition, the contributions by species 
which absorb at the aerosol channel wavelength (but are not 
the target measurement) are also removed from the residual 
extinction particularly if they are measured by HALOE at 
another wavelength. This is also denoted in Table 1 . One im- 
portant exception is that NO2 absorption is not removed from 
the aerosol extinction coefficient measurement at 3.40 pm 
(E. Remsberg, personal communication, 2012). As will be 
demonstrated below, this leaves an artifact in the aerosol ex- 
tinction coefficient values at this wavelength that is notewor- 
thy after 1995. Some species that are not directly measured 
by HALOE (or not measured at all relevant altitudes) are re- 
moved from the residual extinction using climatological val- 
ues for the interfering species. The climatologies used in this 
process were constructed as a part of the Upper Atmosphere 
Pilot Database (UAPD). HALOE makes use of N2O and CH4 
(from near the tropopause and below) from this dataset which 
are based on data from the Stratospheric and Mesospheric 
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Tropics; HALOE 3.46-micron 


TroDics: SAGE II 1020-nm 


Year 


Tropics; HALOE 3.40-micron 


Tropics; HALOE 5. 26-micron 


Fig. 2. This figure shows the 3 -month running median of aerosol 
extinction coefficient between 10° S and 10° N for the 4 HALOE 
channels and the SAGE II 1020-nm measurements. 




Median Relative Reported Uncertainty Median Relative Standard Deviation 


Fig. 3. This figure shows the median relative uncertainty in the four 
HALOE aerosol extinction channels between 1 0° S and 1 0° N and 
1996 and 2005 (left) and the median relative standard based on a 
3 -month running (similar to the analysis shown in Fig. 2) in the 
four HALOE aerosol extinction channels and the SAGE II 1020- 
nm aerosol extinction channel between 1 0° S and 1 0° N and 1 996 
and 2005 (right). 


Sounder (SAMS) mission in 1979 and balloon-based profiles 
(Jackman et al., 1989; Jones and Pyle, 1984). 

Many climatologies of SAD have been derived from 
SAGE II and HALOE independently (Thomason et al., 1997; 
Steele and Turco, 1997; Wang et al., 1989; Yue, 1999; 
Hervig et al., 1998). The goal of this paper is to evaluate the 
HALOE aerosol extinction coefficient data as the first nec- 
essary step in a process that will eventually allow the pro- 
duction of a combined HALOE/SAGE II aerosol climatol- 
ogy of derived aerosol products including surface area den- 
sity (SAD) and a mechanism to produce complete aerosol 
extinction/absorption spectra suitable for use in climate mod- 
eling. There are only very limited directly comparable mea- 
surements to HALOE aerosol extinction coefficient mea- 
surements (e.g., Cryogenic Limb Array Etalon Spectrome- 
ter (CLAES) through May 1993). As a result, past evalua- 
tions of HALOE aerosol extinction coefficient measurements 
have been focused on the strongly volcanic periods prior to 
1995 and measurements from optical particle counters (e.g., 
Hervig et al., 1996) or dependent on using HALOE-derived 
aerosol size distributions to produce comparable bulk prop- 
erties like SAD or extinction at other visible wavelengths 
(e.g., Hervig and Deshler, 2002). Herein, I attempt to eval- 
uate each HALOE aerosol channel with as little dependence 
on assumptions regarding the underlying size distribution as 
possible. I do not consider the HALOE aerosol size distribu- 
tion parameters or bulk parameters derived from those values 
since they are not relevant to the final goal of this study. In 
any case, the HALOE size distribution fits reported in the 
official data files are based on multiple HALOE aerosol ex- 
tinction coefficient values and thus may mask or exacerbate 
issues at individual wavelengths. 


2 General morphology of the HALOE aerosol 
observations 

Figure 2 shows 3 -month running median depictions of 
aerosol extinction coefficient at 2.45 (a), 3.40 (b), 3.46 (c), 
and 5.26 pm (d) for 10° S to 10° N in units of 1 km -1 . Note 
that, while the bulk of the following analysis focuses on the 
tropics, the conclusions HALOE aerosol extinction coeffi- 
cient data quality are independent of latitude. In this case, 
I have interpolated the HALOE extinction to a 0.5 -km grid 
that matches the SAGE II measurement altitudes. Medians 
are used as a crude cloud filter; the median-based results of 
the analysis are essentially identical to mean-based results 
above 1 8 km (where clouds are very rare) and are only mod- 
estly affected by cloud presence well into the troposphere. I 
have also included the SAGE II 1020-nm aerosol extinction 
coefficient analysis analyzed in an identical manner though 
a more formal cloud identification scheme for SAGE II is 
available (Kent et al., 1993). Figure 3 shows a comparison of 
the reported relative measurement uncertainty (Fig. 3a) and 
the relative observed variability (Fig. 3b) for tropical mea- 
surements for (10° S to 10° N) at each of the four HALOE 
aerosol extinction coefficient measurement channels for the 
period 1996 through 2005 where aerosol levels are relatively 
stable. In this case, I show the median value for each param- 
eter coming from the 3 -month depictions shown in Fig. 2. 
In the tropics, low zonal variability is expected during qui- 
escent periods and often observed for a wide collection of 
stratospheric components, particularly above 20 km. As a re- 
sult, the observed variability in this region can be used as 
a rough stand-in for measurement noise or at least an up- 
per bound on these values (e.g., Cunnold et al., 1984). For 
the 3.40, 3.46 and 5.26 pm channels, the median reported 
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measurement uncertainty slowly increases from about 10% 
at 20 km to 20-30 % at 30 km. The observed relative variabil- 
ity in these channels stays roughly in the 1 0 to 20 % range 
through the 20 to 30 km range or comparable to or slightly 
larger than the reported uncertainties. As a result, it appears 
that the reported measurement uncertainties for these chan- 
nels are reasonable estimates of the precision of these mea- 
surements. On the other hand, the reported uncertainty for 
the 2.45 pm channel varies from about 20 % at 20 km to over 
50% at 30 km whereas the observed variability runs from 
15 % at 20 km to 25 % at 30 km. This difference suggests that 
the reported uncertainties overestimate the random compo- 
nent of the error budget (precision) for this channel by about 
a factor of 2. 

While there are significant differences in the aerosol ex- 
tinction kernel (Fig. 1) among the 1020nm, 2.45 pm, and 
the longer wavelength HALOE measurement wavelengths 
(where the kernels are essentially parallel to each other), 
grossly similar behavior in all five analyses should be ob- 
served. In fact, this is broadly what is shown in Fig. 2 (where 
contours are given in logio of extinction). There is an intense 
aerosol layer, centered near 20 km associated with the 1991 
Pinatubo eruption that slowly relaxes toward non-volcanic 
levels by the late 1990s and the aerosol extinction coeffi- 
cient at all wavelengths is fairly stable thereafter. Despite the 
known NO 2 contaminant issue, the 3.40 pm HALOE channel 
is qualitatively very similar to the SAGE II 1020-nm channel 
with a similar structure and duration for the Pinatubo aerosol 
layer, similar seasonal structure in the clean period and ev- 
idence for the impact of two minor eruptions in the 2000 s 
(Ruang in 2002 and Manam in 2005). The 3.40 pm data does 
suggest a more rapid increase in aerosol extinction with de- 
creasing altitude in the upper troposphere (below ~ 17 km) 
than does SAGE II. This may reflect the NO 2 contamination 
or another issue. Since clouds typically exhibit extinction 
coefficients that are many times larger than aerosol, mixed 
views of cloud and aerosol often take on the characteristics 
of cloud events even when the relative fraction of the field 
of view occupied by cloud is relatively small. As a result, the 
more rapid increase in extinction coefficient in at 3.40 pm rel- 
ative to the SAGE II 1020 nm channel may also reflect an in- 
creased spreading of “cloud” observations to higher altitudes 
given the larger field of view for HALOE (1.6 km) relative to 
SAGE II (0.5 km) (e.g., Kent et al., 1997). As shown in Fig. 3, 
SAGE II extinction coefficient relative standard deviation in- 
creases more rapidly below 18 km than do HALOE values. 
This is also consistent with a FOV differences between two 
instruments as relatively infrequent SAGE II cloud observa- 
tions between 1 6 and 1 8 km have a strong impact on the ob- 
served standard deviation but virtually no impact on the me- 
dian extinction coefficient itself. 

The 3.46 pm channel is very similar to the 1020 nm SAGE 
II channel and the 3.40 pm channel and shows a similar 
Pinatubo layer, seasonality in the clear period, and both mi- 
nor volcanic events. In fact, based on the similarity of their 



-25 


%2Q 


1992 



Fig. 4. Median difference between HALOE aerosol extinction coef- 
ficient at 3.40 and 3.46 pm scaled by a factor of 10 5 and in units of 
km -1 . Line contours are drawn at —0.8, —0.6, —0.4, —0.2, 0, 0.2, 
0.4, 0.6, 0.8, 1.0, 3, 5, 7, 9, 20, 40, 60, 80, 100 with an extra color 
contour between each isopleth. 


extinction coefficient kernels, the 3.40 and 3.46 pm channels 
should be virtually identical except for the affect of the afore- 
mentioned NO 2 contamination at 3.40 pm. However, while 
Fig. 2 shows that, as expected, aerosol extinction at 3.46 pm 
is somewhat less than at 3.40 pm in the stratosphere, the dif- 
ference become much larger at lower altitudes as 3.46 pm 
aerosol extinction coefficient decreases rapidly with distance 
below 16 km falling below 10 -5 km -1 (and continuing to de- 
crease) by 12 km (compared to values around 2 x 1(T 4 km" 1 
at 3.40 pm). Figure 4 shows the mean difference between the 
3.40 and 3.46 pm aerosol extinction coefficient data in north- 
ern and southern mid-latitudes and in the tropics. In the mid- 
latitude stratosphere, the differences show a fairly strongly 
annual cycle peaking in the summer and exhibiting a mor- 
phology and magnitude that is roughly consistent with the 
known NO 2 contamination at 3.40 pm. Below 20 km, Fig. 4 
shows that the divergence between the 3.40 and 3.46 pm 
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aerosol extinction coefficient data begins as high as 19 km 
in the tropics and as high as 1 7 km in mid-latitudes or at or 
slightly above the tropopause. There is clearly a defect in 
the 3.46 pm data at lower altitudes and the character of the 
anomaly may suggest the over correction of an interfering 
gas species (overestimating extinction due to absorption) due 
to either a biased concentration (most likely a climatology) 
or spectroscopy for an interfering species in this band. This 
abrupt transition occurs near the altitude at which methane 
clearing transitions from HALOE observations (in the strato- 
sphere) to use of the UAPD methane values (in the tropo- 
sphere). Therefore, UAPD tropospheric methane is a poten- 
tial source of 3.46 pm aerosol extinction coefficient bias. 

Aerosol extinction coefficient at 2.45 and 5.26 pm show 
similar behavior in the stratosphere as the channels at 3.40 
and 3.46 pm. An exception to this is that the minor erup- 
tions in the 2000s are not clearly detectable in the 2.45 pm 
data. At lower altitudes, both channels show a strong in- 
crease with decreasing altitude below 20 km that does not 
conform to what is observed at either 1020nm by SAGE II 
or at 3.40 pm by HALOE. It is not likely due to clouds as 
clouds should influence measurements at all HALOE aerosol 
wavelengths in a similar manner (at least in terms of how it 
is manifested as opposed to actual extinction coefficient val- 
ues). For at least the 5.26 pm channel, there a clue as to the 
source of this anomaly in an annual cycle between 15 and 
23 km whose phase and tilt in time that is very analogous to 
the water vapor tape recorder (Mote et al., 1996). Taken with 
the rapid increase at and below the tropopause, the behavior 
of the 5.26 pm channels data suggests a water vapor-based 
artifact in this channel though the artifact could be derived 
from another gas species that follows a temporal/spatial pat- 
tern similar to water vapor. HALOE ’s 2.45 pm aerosol ex- 
tinction coefficient data has a similar rapid increase below 
15 km as shown by the 5.26 pm channel. However, there is 
no evidence of a tape recorder like structure in the strato- 
sphere. There is some weak absorption by ozone and N 2 O 
within the 2.45 pm channel filter coverage and it is not clear 
that either is cleared using measured ozone or UAPD N 2 O. 


3 Anomaly analysis 

Quantifying the overall quality of the HALOE aerosol ex- 
tinction coefficient measurements, particularly in the rela- 
tively clean period after 1996, is limited by the lack of di- 
rectly comparable measurements and has, in the past, been 
based on the conversion of in situ measurements by the Uni- 
versity of Wyoming Optical Particle Counter (OPC) (Hervig 
et al., 1993) and by conversion of the HALOE measure- 
ment themselves into extinction at other wavelengths or other 
derived products (e.g., SAD) using the size distributions 
derived using multiple HALOE measurement wavelengths 
(Hervig et al., 1995, 1998; Hervig and Deshler, 2002). Since 
the goal of this work is to assess individual HALOE chan- 


nels at near background levels, derived quantities based on 
size distribution fits are not well suited to this evaluation. 
Also, while a single mode log-normal is often used in mod- 
eling the stratospheric aerosol size distribution and is a rea- 
sonable approximation for observed size distributions, multi- 
mode and more complex aerosol size distributions are com- 
monly observed (e.g., SPARC, 2006) and observed aerosol 
extinction coefficient are often sensitive to the details of the 
size distribution. While it is not totally possible to eliminate 
some sort of a priori model for the size distribution, I have 
attempted to minimize the sensitivity of the following evalu- 
ation to those assumptions and make no attempt to infer the 
underlying size distribution. Instead, I use a family of single 
mode log-normal size distributions (ranges in the mode ra- 
dius and width) to roughly carve out the space in which ratios 
of pairs of different wavelengths of aerosol extinction coeffi- 
cient should exist. Using ratios eliminates sensitivity to total 
aerosol number density and the size distribution can be con- 
sidered a smoothing function for the aerosol extinction ker- 
nels. As a result, “ratio space” defined through log-normal 
functions should be fairly representative of any reasonable 
stratospheric aerosol size distribution. I define a good match 
to be whenever the observed ratio pairs lie within or near 
the model space (within ~ 20 % of the bounds defined by 
the models). I do not require good matches from one chan- 
nel to lie near the equivalent location within the log-normal 
space (i.e., the same implied log-normal size distribution) in 
a further attempt to minimize the impact of aerosol size dis- 
tribution assumptions. 

Figure 5 shows the distribution of the 3 -month running 
medians (shown as the symbol “+”) for the ratio of a 
HALOE aerosol extinction coefficients relative to the SAGE 
II 1020 nm value versus the SAGE II 525 to 1020-nm aerosol 
extinction coefficient ratio (derived using the methodology 
used to produce Fig. 2. Also shown are log-normal models 
results shown as lines for constant size distribution width at 
values of 1.2, 1.4, 1.6, and 1.8. The SAGE II 525 to 1020 nm 
aerosol extinction coefficient ratio can be more or less trans- 
lated as time since the record begins in 1991 when aerosol 
showed a small extinction ratio (implying large aerosol) and 
which increases (~ monotonically) into the 2000s. This fig- 
ure shows that the early data at this altitude for all four 
HALOE channels are in fairly good agreement with the 
SAGE II ratios when the SAGE II ratio is less than 2. This 
occurs early in the HALOE record when extinction is dom- 
inated by the Pinatubo event and are the largest observed 
during the HALOE mission lifetime. Following that period, 
the ratios using 3.40 and 3.46 pm HALOE channels remain 
within or close to the single mode-lognormal space through- 
out their records. On the other hand, both the 2.45 and 
5.26 pm ratios consistently increase to where they lay as 
much as a factor of 2 (at 5.26 pm) or 3 (at 2.45 pm) outside 
the area bound by the single-mode lognormal model. While 
some leeway must be granted for the potential for deficien- 
cies in the SAGE II data and the existence of more complex 
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Fig. 5. These figures show the HALOE 2.45 (a), 3.40 (b), 3.46 (c) 
and 5.26 pm (d) to SAGE II 1020-nm aerosol extinction coefficient 
ratio relative to the SAGE II 525 to 1020-nm aerosol extinction ratio 
(shown as “+” symbols) computed following the 3 -month running 
medians as shown in Fig. 2 for the Tropics and at 20 km. Data from 
prior to 1997 is shown in blue. The lines represent model relation- 
ships for sulfate aerosol at stratospheric temperatures based on a 
single mode log-normal with varying mode radius (along the lines) 
and width (from 1.2 for the solid line to 1.8 for the dash-dot line). 


aerosol size distributions, it is virtually impossible to con- 
struct any size distribution that would match the observed be- 
havior at 2.45 and 5.26 pm relative to SAGE II observations 
particularly without simultaneously breaking the positive be- 
havior observed at the 3.40 and 3.46 pm channels. This figure 
implies that the 2.45 and 5.26 pm HALOE aerosol extinction 
coefficient measurements are not consistent with the HALOE 
3.40 and 3.46 pm values as well as the SAGE II measure- 
ments. 

Selecting the 3.46 pm channel as the nominal standard for 
HALOE aerosol extinction coefficient channels (in much the 
same way as SAGE II’s 1020nm channel is treated as the 
paramount member of the SAGE II ensemble), Fig. 6 shows 
an analysis similar to that in Fig. 5 using the 3.46 pm to 
1020nm ratio channel as the independent variable relative 
to ratios of 2.45, 3.40 and 5.26 pm relative to 3.46 pm as the 
dependent variable. This reduces the dependence on com- 
parisons to SAGE II measurements particularly for channels 
at 3.40 and 5.26 pm where the model results show virtually 
no variation with the 3.46 pm to 1020 nm ratio. On the other 
hand, the 2.45 to 3.46 pm extinction coefficient ratio remains 
sensitive to the use of the SAGE II data since sulfate aerosol 
is only weakly absorbing at 2.45 pm and the resulting ex- 
tinction kernel (shown in Fig. 1) is a mix of features com- 
mon to SAGE II aerosol kernels and the longer wavelength 
HALOE channels where sulfate absorption is very strong. 
As in the previous figure the 3.46 pm to 1020nm extinc- 
tion coefficient ratio is also a rough stand in for increasing 



Fig. 6. These figures show the HALOE 2.45 (a), 3.40 (b) and 
5.26 pm (c) to 3.46 pm aerosol extinction ratio relative to the 
HALOE 3.46 pm to SAGE II 1020-nm aerosol extinction ratio 
(shown as “+” symbols) computed following the 3 -month running 
medians as shown in Fig. 2 for the Tropics and at 20 km. Data from 
prior to 1997 is shown in blue. The lines represent model relation- 
ships for sulfate aerosol at stratospheric temperatures based on a 
single mode log-normal with varying mode radius (along the lines) 
and width (from 1.2 for the solid line to 1.8 for the dash-dot line). 


time and decreasing absolute extinction coefficient. In Fig. 6, 
the offset of the 3.40 pm channel from its expected relation- 
ship with the 3.46 pm channel (and virtually independent of 
the SAGE II 1 020-nm extinction coefficient) averages about 
+10% while the 5.26 to 3.46 pm (expected to be similarly 
independent of the 1 020 nm extinction) varies from +20 % 
at smallest ratios (i.e., the highest aerosol loading) to more 
than 100 % at larger ratios. Both of these figures yield inter- 
pretations consistent with those drawn from Fig. 5. On the 
other hand, the 2.45 to 3.46 pm figure suggests that the data 
is more consistent with the 3.46 pm and 1020-nm extinction 
coefficient ratio than would be inferred from the more SAGE 
II-based comparison shown in Fig. 5 for at least the first half 
of the data set before substantial departures from the model 
predictions appear. However, the relationship among the data 
shows almost no dependence on 3.46 pm to 1020nm extinc- 
tion coefficient ratio unlike the strong dependence based on 
the model results. As with the 5.26-pm channel, results from 
the 2.45- pm channel appear incorrect except at the highest 
levels of aerosol extinction coefficients observed. 

I have computed an aerosol extinction coefficient anomaly 
as a function of altitude using the single mode log-normal 
model with a width of 1.6 (the dashed lines in Figs. 5 
and 6) for all four HALOE aerosol channels relative to 
the SAGE II ratios shown in Fig. 5 and for the 2.45, 3.40 
and 5.26 pm channels relative to the HALOE 3.46 pm and 
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Fig. 7. The figure shows the inferred aerosol extinction coefficient anomalies (the difference between the single mode log-normal size 
distribution model and the observed channel values) based on the SAGE II 525 to 1020-nm model shown in Fig. 5 for 20 km (a) and based 
on the HALOE 3.46 pm to SAGE II 1020-nm model (b and c) shown in Fig. 6 for a size distribution width of 1.6. The inset in (a) shows the 
1020-nm anomaly below 17 km with an expanded extinction coefficient scale (in km -1 ). 


SAGE II 1020-nm channel ratios shown in Fig. 6. In Fig. 5, 
the anomaly is computed as the SAGE II 1020-nm extinc- 
tion coefficient times the difference between the model pre- 
diction for the HALOE extinction coefficient values (rela- 
tive to SAGE II 1020-nm extinction and based on the ob- 
served SAGE II 525 to 1020-nm extinction) and the observed 
HALOE extinction relative to the SAGE II 1020-nm extinc- 
tion coefficient. For this analysis I have limited the altitude 
range to 15 to 30 km and show the median differences for 
1996 through the end of the record and latitudes from 10° S 
to 10° N. Given the limitations of the anomaly calculation, 
there is some absolute uncertainty that may vary as a func- 
tion of time and altitude (as the true underlying size distribu- 
tion changes) so that the differences shown in Fig. 7 should 
be considered only estimates of the potential for bias in any 
of the channels. Also, since I have based this analysis on 
the 3.46 pm channel care must be used in interpreting the 
analysis below 19 km where there is obvious degradation in 
this channels performance. Despite this, I used the 3.46 pm 
channel as the base HALOE extinction measurement because 
it lacks the known NO 2 -based bias of the 3.40 pm channel 
throughout the depth of the profile and other channels are 
far more suspect. Figure 7a shows the anomalies in the all 
four channels relative to the SAGE II-based predictions and 
demonstrates a consistent picture between 20 and 30 km of 
small anomalies (~ zb 10 -6 km -1 ) for 3.40 and 3.46 pm but 
significantly larger ones at 2.45 and 5.26 pm (~ 10 -5 km -1 ). 
However, all four channels show much larger deviations from 
the model predictions below 1 5 km (shown in the inset on 
Fig. 7a) and with values larger than 10 -3 km -1 at 5.26 pm. 
Figure 7b and c shows the anomalies based on the 3.46 pm 
to 1020-nm aerosol extinction ratio in absolute and relative 


sense. The differences between the 3.40 and 3.46 pm extinc- 
tion coefficients are on the order of 10 -6 km -1 or about 10 % 
of the overall measurement between 20 and 30 km. This dif- 
ference is fairly consistent with the expected NO 2 absorption 
at 3.40 pm and suggests that these two channels would be 
in excellent agreement in this altitude range if the effects of 
NO 2 absorption were removed from the 3.40 pm channel. On 
the other hand, the channels at 2.45 and 5.26 pm show differ- 
ences on the order of 0.5 to 1 x 10 -5 km -1 which represent a 
20 to 50 % anomaly at 5.26 pm and 60 to 95 % at 2.45 pm. It 
is possible that more complicated size distributions and more 
diverse composition (e.g., ice or organics) would account for 
some of the apparent deficiencies observed herein, however, 
I believe that the bulk of the problems are HALOE data qual- 
ity issues. 

It is beyond the scope of this paper to fully diagnose and 
repair deficiencies in HALOE aerosol extinction coefficient 
data. Nonetheless, given inferences I have already drawn re- 
garding the potential for the residuals to arise from the failure 
to effectively remove interfering species from the aerosol re- 
trievals, I have made further efforts to identify potential can- 
didates as sources for the observed data quality issues. To do 
this, I have computed the correlation and its uncertainty be- 
tween the estimated aerosol extinction coefficient anomalies 
and the seven gas species reported by HALOE. As in Fig. 7, 
I have limited the analysis to 15 to 30 km and for HALOE 
events from 1996 to the end of the record and latitudes from 
10° S to 10° N. It should be kept in mind that correlation be- 
tween the estimated aerosol anomalies and any gas species 
may imply a causal relationship between these parameters or 
a mutual correlation to another, perhaps unmeasured, param- 
eter. Figure 8 shows the correlation of the 3.46 pm anomaly 
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Fig. 8. Correlation of 3.46 jam anomalies based on the SAGE II 525/1020 model with HALOE measured species. 
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Fig. 9. Correlation of 2.45 pm anomalies based on the HALOE 3.46 pm model with HALOE measured species. 
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Fig. 10. Correlation of 3.40 jam anomalies based on the HALOE 3.46 jam model with HALOE measured species. 


(computed using the SAGE II ratios as in Fig. 5) and the 
seven HALOE gas species. In this case, above 20 km, no con- 
sistent evidence for correlations between any gas species and 
the anomaly is observed though some small but significant 
correlations with ozone, CH4, and HF can be noted. Below 
20 km, some moderate and significant negative correlations 
with HC1, ozone, and H2O are also seen near 17 km though 
they have mostly returned to near zero at 1 5 km and so it is 
not clear if they have any relationship to the bias noted in the 
3.46 pm channel at and below this altitude. 

Figure 9 shows the correlation of the 2.45 pm anomaly 
relative to the HALOE gas species measurements. In this 
case, we see a clear positive correlation between the anomaly 
and HF, which is measured at 2.45 pm, between 18 and 
28 km which suggests an incomplete removal of HF from 
the aerosol residual. In addition, significant correlations are 
noted with ozone between 20 and 25 km (where ozone is a 
maximum), and NO and NO2 between 23 and 28 and be- 
low 20 km (with a change in sign). Ozone weakly absorbs at 
this wavelength and is not listed among the species cleared in 
the retrieval process (Hervig et ah, 1996) and so the aerosol 
anomalies may plausibly be related to the failure to clear 
this species. On the other hand, NO and NO2 do not absorb 
at this wavelength so the significant correlation may be re- 
lated to another species or some other process that is corre- 
lated with NO and/or NO2. Figure 10 shows the correlation 
of the 3.46 pm anomaly with the HALOE gas species mea- 
surements. In this case, NO2 is positively correlated with the 


aerosol anomaly throughout the 1 5 to 20 km range and re- 
inforces the impact of the failure to clear NO2 at this wave- 
length. NO is similarly correlated but this is most likely due 
to the correlated between NO and NO2 themselves. Ozone 
has a similar correlation, but opposite in sign, which may 
again be the result of correlation of ozone and NO2 rather 
than a signature of residual ozone in the aerosol data though 
that cannot be totally ruled out. Water vapor also shows some 
significant positive correlations with some fairly complex 
vertical structure. Since water is removed using HALOE ob- 
servations, it is not clear why this correlation exists. Finally, 
there are some surprisingly strong correlations with HF (in 
particular) and HC1 although neither has appreciable absorp- 
tion at this wavelength. At this time, I do not fully understand 
the source of this correlation. Figure 1 1 shows the correlation 
of the 5.26 pm anomaly with the HALOE gas species mea- 
surements. Here, I find strong correlation with water vapor 
below 23 km which seems consistent with the observation of 
a “tape recorder” like feature in the aerosol extinction coeffi- 
cient data at this wavelength. I also found strong correlations 
with HF, NO and NO2 above 23 km although none of them 
absorb at 5.26 pm. It is possible that the strong correlation 
with NO and NO2 is related to either the failure to remove or 
an ineffective removal of N2O which has significant absorp- 
tion at 5.26 pm. Like correlations at 2.45 pm, I do not fully 
understand the source of the correlations exhibited by this 
channel. 
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Fig. 11. Correlation of 5.26 pm anomalies based on the HALOE 3.46 pm model with HALOE measured species. 


4 Application of HALOE data 

The evaluation of HALOE aerosol extinction coefficient data 
found that one at 3. 46- pm channel, is of sufficient qual- 
ity (without further corrections) to use in an evaluation of 
aerosol properties (> 19 km). While multiple channels would 
be preferred, in combination with multiple SAGE II chan- 
nels, the single channel adds significantly to the information 
contained in the ensemble and holds the potential to sub- 
stantially improve estimates of surface area density (SAD) 
and other bulk properties. Thomason et al. (2008) showed 
that a notional measurement of aerosol volume substantially 
reduced the uncertainty inherent in visible wavelength-only 
measurements as well as modestly increasing the overall es- 
timates for SAD. Following the method used in that paper, 
I have produced a simple model to combine SAGE II ex- 
tinction coefficient measurements at 525 and 1020nm with 
those at 3.46 pm from HALOE. In this approach, I use the 
525 and 1020nm extinction coefficient values to estimate 
a monomodal “large” mode aerosol that reproduces the ex- 
tinction at both wavelengths exactly. SAD calculated from 
this mode is the minimum possible SAD that is consistent 
with those measurements (Thomason et al., 2008). Gener- 
ally, the radius of this mode is between 0.2 and 0.3 pm and 
the number densities associated with this fit are on the order 
of 1 cm -3 and much less than the nominal stratospheric value 
of ~ 10 cm -3 . The extinction at 3.46 pm implied by this fit 
is also much less than observed by HALOE and effectively 


demonstrates the insensitivity of visible wavelength aerosol 
extinction to small aerosol (< 0.1 pm). I then introduce a sec- 
ond monomode which has a number density which brings the 
total number density to either 5 or 15 cm -3 . The size of this 
mode is selected such that the extinction at 3.46 pm is repro- 
duced but do not require either SAGE II extinction to remain 
fixed. The radii of this mode are usually less than 0. 1 pm. I 
find the extinction computed using both modes reproduces 
the 1020 nm extinction within a few percent of the measured 
value. The impact at 525 nm is somewhat larger and is usu- 
ally between 1 and 2 times the 525-nm measurement un- 
certainty (10 to 20%). The impact at 525 nm could be re- 
duced with a slightly more sophisticated model. The average 
of the two fits is used as the “final” SAD product with the 
profiles based on the 5 and 15 cm -3 representing the lower 
(from 15 cm -3 ) and upper (from 5 cm -3 ) bounds for SAD. 
Figure 12 shows the results for the tropics in March 1999 
using the method of Thomason et al. (2008) (based SAGE 
II alone) with an estimated range of close to a factor of 3 
throughout the profile. The figure also shows results using the 
method described above. These results lie within the SAGE 
II-only result range but generally the mean value is about 
30 % higher. Most significantly the range of SAD values is 
reduced by more than a factor of 3. Obviously, this model 
does not (and is not intended to) produce physically realistic 
aerosol size distributions but is instead intended to produce a 
viable range for retrieved bulk properties like SAD consistent 
with observed extinction values but minimally dependent on 
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Fig. 12. Comparison of SAD estimates for March 1999 in the trop- 
ics (10° S to 10°N) based on SAGE II-only (black) based on the 
method of Thomason et al. (2008) and using SAGE II 525 and 1020- 
nm extinction coefficient data along with that from the HALOE 
3.46 pm based on a modified form of the Thomason et al. (2008) 
method (red). 

assumptions regarding the underlying size distribution. In 
that regard, this result demonstrates the potential positive im- 
pact of employing a mixed visible/infrared aerosol extinction 
ensemble in the derivation of climate/chemically important 
aerosol parameters. 


5 Conclusions and recommendations 

In summary, based on the analyses above for the time pe- 
riod after 1996, 1 have concluded that HALOE’s 3.46 pm is 
of good quality above 19 km and suitable for scientific ap- 
plications above that altitude. However, I find that it is of 
increasingly suspect quality at lower altitudes and should not 
be used below 17 km under any circumstances. The 3.40 pm 
is biased by about 1 0 % throughout the stratosphere due to 
the failure to clear NO 2 but otherwise appears to be a high 
quality product down to 15 km. Its behavior below 15 km is 
more cause for concern and should be used with caution. Fi- 
nally, I find that both the 2.45 and 5.26 pm aerosol extinction 
coefficient measurements are clearly biased and should not 
be used for scientific applications after the strongest parts of 
the Pinatubo period. It is clear that the 3.40 pm aerosol ex- 
tinction coefficient measurements can be improved through 
the inclusion of an NO 2 correction and could, in fact, end up 
as the highest quality overall HALOE aerosol extinction co- 
efficient measurement. This could be done using HALOE’s 
NO 2 product in the production of a future version or done 
using a relatively simple model for NO 2 absorption (e.g., 
MODTRAN) as a post- version correction. On the other hand, 


fully understanding the deficiencies of the 2.45 and 5.26 pm 
channels probably requires access to the HALOE process- 
ing code and supporting data sets and more information re- 
garding these than I have been able to acquire at this time. 
Nonetheless, careful examination of the UAPD would seem 
like an important activity. The unique characteristics of the 
2.45 pm aerosol kernel (size sensitivity with some absorp- 
tion) make it potentially very valuable in inferring the un- 
derlying characteristics of the aerosol size distribution. Im- 
provements to this channel should be considered a high pri- 
ority. On the other hand, since 5.26 pm aerosol extinction co- 
efficient measurements do not offer information that is par- 
ticularly unique relative to other measurements, future im- 
provements to this channel do not have the same value. Fi- 
nally, a simple model to demonstrate the power of mixed vis- 
ible/infrared aerosol extinction coefficient ensembles for the 
retrieval of bulk aerosol properties demonstrates that a com- 
bined HALOE/SAGE II aerosol climatology is feasible and 
may represent a substantial improvement over independently 
derived data sets. 
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